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Abstract 

The explicit semiclassical treatment of the logarithmic perturbation theory for 
the bound-state problem for the spherical anharmonic oscillator and the screened 
Coulomb potential is developed. Based upon the ^-expansions and suitable quanti- 
zation conditions a new procedure for deriving perturbation expansions is offered. 
Avoiding disadvantages of the standard approach, new handy recursion formulae 
with the same simple form both for ground and excited states have been obtained. 
As examples, the perturbation expansions for the energy eigenvalues of the quartic 
anharmonic oscillator and the Debye potential are considered. 

1 Introduction 

The main task in application of the quantum mechanics is to solve the Schrodinger equa- 
tions with different potentials. Unfortunately, realistic physical problems can practically 
never be solved exactly. Then one has to resort to some approximations. Most widely 
used among them is the perturbation theory. However, the explicit calculation with 
the Rayleigh - Schrodinger perturbation theory, described in most quantum mechanics 
textbooks, runs into the difficulty of the summation over all intermediate unperturbed 
eigenstates. To avoid this difficulty, various alternative perturbation procedures have been 
proposed [1-13]. 

Nevertheless up to now, one of the principal approximation techniques is the loga- 
rithmic perturbation theory [14-21]. Within the framework of this approach, the con- 
ventional way to solve a quantum-mechanical bound-state problem consists in changing 
from the wave function to its logarithmic derivative and converting the time-independent 
Schrodinger equation into the nonlinear Riccati equation. 

In the case of ground states, the consequent expansion in a small parameter leads to 
handy recursion relations that permit us to derive easily the corrections to the energy as 
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well as to the wave function for each order. However, when radially excited states are con- 
sidered, the standard technique of the logarithmic perturbation theory becomes extremely 
cumbersome and, practically, inapplicable for describing higher orders of expansions. 

Recently, a new procedure based on specific quantization conditions has been proposed 
to get series of the logarithmic perturbation theory via the ^-expansion technique within 
the framework of the one-dimensional Schrodinger equation [22] . Avoiding the disadvan- 
tages of the standard approach, this straightforward semiclassical procedure results in 
new handy recursion formulae with the same simple form both for the ground state and 
excited states. 

The object of the present work is to extend the above mentioned formalism to the 
bound-state problems within the framework of the three-dimensional Schrodinger equa- 
tion with central potentials, such as the anharmonic scillator potential and the screened 
Coulomb one, which are widely used in practice. 

2 Basic concepts of the method 

We study the bound-state problem for a non-relativistic particle moving in a central 
potential admitted bounded eigenfunctions and having in consequence a discrete energy 
spectrum. Let us therefore consider the reduced radial part of the Schrodinger equation 

with the effective potential having only one simple minimum. 

Following usual practice, we apply the substitution, C(r) = hU'(r)/U(r), accepted in 
the logarithmic perturbation theory and go over from the Schrodinger equation ([1]) to the 
Riccati equation 

W{r) + C 2 {r) = - ^ + 2mV{r) - 2mE. (2) 

According to our assumption, we are seeking the eigenvalues and the eigenfunctions 
of this equation explicitly in a semiclassical manner with series expansions in powers of 
the Planck constant 

oo oo 

E=Y, C(r) = ]T C k (r)h k , (3) 

where the order in h of these quantities, i.e. the values of kl and k2, should be defined 
preliminary. 

As in the standard approach, the substitution of these expansions into the Riccati 
equation leads to the simple recursion system. This system can be solved successively in 
the case of ground states, while the description of the excited states has some problems 
with taking into account the nodes of wave functions. For avoiding these problems, 
we shall attempt to use the quantization condition and the formalism of the theory of 
functions of a complex variable. 

Remind that, in the complex plane, a number of zeros iV of a regular function inside a 
closed contour is defined by the principle of argument known from the complex analysis. 
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Being applied to the logarithmic derivative, C(r), it means that 



— J>C(r)dr = — YV J)C k (r)dr = hN. (4) 



2m J 2m 

J fc=0 

This quantization condition is exact and is widely used for deriving higher-order cor- 
rections to the WKB-approximation [23,24] and the 1/iV-expansions [25-27]. There is, 
however, one important point to note. Because the radial and orbital quantum numbers, 
n and I, correspondingly, are specific quantum notions, the quantization condition (J3J) 
must be supplemented with a rule of achieving a classical limit for these quantities. It is 
this rule that stipulates the kind of the semiclassical approximation. 

In particular, within the framework of the WKB-approach the passage to the classical 
limit is implemented using the rule 

H — > 0, n — > oo, I — > oo, hn — const, hi = const, (5) 

whereas the 1/iV-expansion requires the condition [25-27] 

h — > 0, n — const, I — > oo, hn — > 0, hi — const. (6) 

The proposed semiclassical treatment of the logarithmic perturbation theory involves 
the alternative possibility: 

h — > 0, n — const, I = const, hn — > 0, hi — > 0. (7) 

With the last rule, the right-hand side of the equation (4) has the first order in h and 
the quantization condition now takes the simple form 

(fc 1 (r)dr = N, -L (fc k (r)dr = 0, k^l. (8) 
2m J 2m J 

However, this definition of the quantization condition is incomplete since we have not 
pointed out the path of integration. We shall now show that the suitable choice of the 
contour of integration and the consequent integration with application of the Cauchy 
residue theorem easily solves the problem of describing radially excited states. 



3 The anharmonic oscillator 

Quantization conditions. The discussion of details of the proposed technique we begin with 
the case of the anharmonic oscillator potential which is given by a symmetric function 
V(r) that can be written as a Taylor series expansion 

V(r) = -muo 2 r 2 + ^ f iT 2i+ \ (9) 

i>l 

In the first place, let us consider the rule (7) of achieving the classical limit from the 
physical point of view. Since hi — » as h — > 0, the centrifugal term, h 2 l (/ + 1) /r 2 , has 
the second order in h and disappears in the classical limit that corresponds to falling a 
particle into the center. This means that a particle drops into the bottom of the potential 
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well as H — > and its classical energy becomes E = mmV(r) = 0. Hence, the series 
expansions in powers of the Planck constant for the energy eigenvalues and the C (r) must 
now read as E = J2T=i ^^h k anc ^ ^( r ) = ^2T=o Ck(j)Ti k . 

Upon inserting these expansions into the Riccati equation ([2]) and collecting coefficients 
of equal powers of h, we obtain the following hierarchy of equations 

C 2 (r) 

C' (r) + 2C Q {r)C 1 (r) 
C[{r) + 2Cv{r)C 2 {r) + Cl{r) 

k 

Cfc-iW +J2 C ^ C ^ r ) = ~ 2mE ^ k>2. (10) 

i=0 

In the case of ground states, this recurrence system can be solved as straightforwardly 
as in the standard approach. For excited states, we intend to take into account the nodes 
of the wave function with the quantization condition (8) for which we must define the 
contour of integration. 

It should be stressed that our technique is quite distinguished from the WKB method 
not only in the rule of achieving a classical limit but also in the choice of a contour of 
integration in the complex plane. With a view to elucidate the last difference let us 
now sketch out the WKB treatment of this bound-state problem. In the complex plane, 
because the potential is described by the symmetric function ([9]), there are two pairs of 
turning points, i.e. zeros of the classical momentum, on the real axis. Therefore we have 
two cuts between these points: in the region r > as well as in the region r < 0. In spite 
of only one cut lies in the physical region r > 0, the contour of integration in the WKB 
quantization condition has to encircle both cuts for the correct result for the harmonic 
oscillator to be obtained [28]. 

In our approach, when a particle is dropping into the bottom of the potential well 
these four turning points are drawing nearer and, at last, are joining together at the 
origin. Consequently, all zeros of the wave function are now removed from both positive 
and negative sides of the real axis into the origin and our contour of integration must 
enclose only this point and no other singularities. 

Further, let us count the multiplicity of a zero formed in the wave function at r = 0. 
Imposed by the requirement of the regularity, the behavior r l+1 as r — > brings the value 
I + 1. The number of nodes in the physical region r > is equal to the radial quantum 
number n. But, because the potential is a symmetric function, the same number of 
zeros must be in the region r < 0, too. Then the total number of zeros inside the contour 
becomes equal to N — 2n + I + 1. 

For evaluation of the contour integrals in the quantization condition (8), let us consider 
the system fflUj) and investigate the behavior of the functions Cfc(r). From the first 
equation of this system, it follows instantly that the Co(r) can be written as 

/ oo \V2 

C (r) = -\2mV(r)] 1/2 = -mujr 1 + X"f ir 2i ) = rVc°r 2i , (11) 

\ i=l / i=0 



2 mV{r) 
-2 mE x 
1(1 + 1) 



- 2 mE 2 , 
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where the minus sign is chosen from the boundary conditions and coefficients C 4 are 
defined by parameters of the potential through the relations 



i-i 

C° = -mw, C° = ^- ( C p C Ip - 2m/, ) , i > 1. (12) 



At the origin, on account of the equality Co(0) = 0, a simple pole arises for the function 
Ci(r), while Ck{r) has a pole of the order (2k — 1). Thus Ck(r) can be represented by 
the Laurent series 

oo 

C k {r) = r l - 2k E Cy\ k > 1. (13) 

j=0 

Finally, with applying the residue theorem, the quantization condition (8) expressed 
explicitly in terms of the coefficients Cf takes the especially simple form 

CU = N6 lik , (14) 

where N = 2n + I + 1 . 

It is this quantization condition that makes possible the common consideration of the 
ground and excited states and permits us to derive the simple recursion formulae. 

Recursion formulae and the example of application. The substitution of the series 
(I12p and ( |T3l) into the system (ITOl) in the case i ^ k — 1 yields the recursion relation for 
obtaining the Laurent-series coefficients of the logarithmic derivative of the wave function 



2C ° 



k— 1 i i 

(s-2k + 2i)c*- 1 + E E ( V' S : I + 2 E c l c U - ^ + l)Mo,. 

(15) 

If z = — 1, by equating the expression (JT5l) for to the quantization condition 
( HM we arrive at the recursion formulae for the energy eigenvalues 

2mE k = -C k k z\ -EE C i C tlp ■ ( 16 ) 

j=0 p=0 

Derived in this way, first corrections to the energy eigenvalues of the spherical anhar- 
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monic oscillator take the form 

1 + 2N (3-2L + 67/) j\ 

tjX = o U ; ' 'a 2 — 2 ' 

2 Am z uj z 

E * = ^4 ((~ 21 + 9 L - 17 V) h 2 + m (15 - 6 L + 10 r?) ^ / 2 ) 

8 m 4 a; 

£ 4 = \ s ( (333 + llL 2 -3L (67 + 86t/) + 3t/ (347+ 125 77)) f\ 
lo m° cu 8 \ 

-6m (60 + 3 (-13 + L) L + 175 77 — 42 L 77 + 55 7/ 2 ) w 2 /1 / 2 
+m 2 (6L 2 - 12L (6 + 57/) + 35 (3 + 2t/ (4 + 7/))) cu 4 , 

E 5 = l + 2 g iV 11 ( (30885 + 909 L 2 — 27 L (613 + 330 77) + 7/ (49927 + 10689 7/)) / x 

128 77l 8 LU LL \ 

-4 m (11220 + 393 L 2 - 6 L (1011 + 475 77) + 77 (16342 + 3129 7/)) u 2 h 2 f 2 
+ 16 m 2 (33 L 2 — L (501 + 190 77) + 63 (15 + 7/(19 + 3 77))) co A h f 3 
+2 m 2 (3495 + 138 L 2 + 4538 77 + 786 rj 2 - 30 L (63 + 26 7/)) c/ / 2 2 



-4m 3 (30L 2 -20L (24 + 7t/) + 63 (15 + 27/ (8 + 77))) u 6 /A 



(17) 



where AT = 2 n + Z + 1, 77 = N (N + 1), L = 1(1 + 1). 

As it was expected, the obtained expansion is indeed the series of the logarithmic 
perturbation theory in powers of the Taylor-series coefficients for the potential function, 
with the first approximation being equal to the energy of the three-dimensional harmonic 
oscillator 

Ex = ( 2n + l + - J u. (18) 



2, 

Thus, the problem of obtaining the energy eigenvalues and eigenfunctions for the 
bound-state problem for the anharmonic oscillator can be considered solved. The equa- 
tions (|T5|) - (|T6|) have the same simple form both for the ground and excited states and 
define a useful procedure for the successive calculation of higher orders of expansions of 
the logarithmic perturbation theory. 

As an example, we examine eigenenergies for the anharmonic oscillator with the po- 
tential 

V(r) = mu 2 r 2 /2 + Ar 4 , A > 0. (19) 
Then the equations (j!7p are rewritten as 

2 / 4m z uJ z 



E _ -(1 + 2AQ (21-9L + 177/) >2 



<m A u! 5 



(333 + 11 L 2 - 3 L (67 + 86 rj) + 3 7/ (347 + 125 77)) 
16 mP oo a 

-(1 + 2 AT) (30885 + 909 L 2 - 27 L (613 + 330 7/) + r/ (49927 + 10689 7/)) 4 
5 ~ 128m 8 cu n ' 

We recall that here N = 2 n + I + 1, rj = N (N + 1), L = 1(1 + 1). 



It is readily seen that the use of the ^-expansion technique does lead to the explicit 
perturbation expansion in powers of the small parameter A. 

In the case of ground states, obtained expansions for the energy eigenvalues coincide 
with those listed in Ref. [29]. In the case of excited states, our corrections coincide with 
corrections up to the second order which are just only calculated in Ref. [30]. 

As it is known, the expansions for the anharmonic oscillator are asymptotic and di- 
verge for any finite value of the parameter A that requires the use of some procedures of 
improving the convergence (for references see [31]). It should be noted, that the proposed 
technique is easily adapted to apply any scheme of the series renormalization [32]. 



4 The screened Coulomb potential 

Quantization conditions. Now let us consider the case of the screened Coulomb potential 
which in common practice has a form 



V(r) = -F(k,t). 



(21) 



where k is a small parameter. 

In what following, we do not single out explicitly the screening parameter, but incor- 
porate it into coefficients V% of the Taylor series expansion of this potential 



V(r) = -J2 



(22) 



i=0 



Note, that after performing the scale transformation r — > h 2 r powers of the screening 
parameter appear in common with powers of Planck's constant squared. Hence, the 
perturbation series must be, as a matter of fact, not only ^-expansions but also the 
semiclassical h 2 -expansions, too. 

In the classical limit, when a particle falls into the center, its energy eventually ap- 
proaches infinity. Hence, the expansions ([3]) must be represented as E = h~ 2 ^2^ =0 Ekh 2k 
and C(r) = h^ 1 ^2^ =0 Ck(r)h k that results in the recurrent system 



C 2 (r) = -2mE , 
C (r)d(r) =m\v(r) 



Ex 



C' 1 (r) + 2C (r)C 2 (r) + C 2 (r) 



1(1 + 1) 



-2mE 2 



-2mE k , k > 2. 



(23) 



which changes only in the first two equations in comparison with (10). 

Now, let us consider the choice of the contour of integration in the quantization rela- 
tion. Since in the classical limit a particle falls into center, the classical turning points 
again draw to the origin and the nodes of the wave function are joined together at r = 0. 
Then, as well as in the anharmonic oscillator case, the contour of integration must enclose 
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only the origin. However, now the nodes of the function come into the origin only from 
the positive region of the real axis. Thus, with the number of nodes and the value I + 1 
included, the total number of zeros in the quantization condition (jlj) becomes equal to 
N = n + l + l. 

Further, from (1231) it appears that Co(r) is the constant and its Taylor-series coefficients 

are 

C° = -v/-2m£ , Cl = 0, (24) 

Owing to the Coulomb behavior of the potential at the origin, the C*i(r) has a simple pole 
at this point, while the function C k (r) has a pole of the order k and may be represented 

as 

oo 

C k (r)=r- k J2cy, (25) 



i=0 



that leads to the known quantization condition 



(26) 



where N = n + I + 1 . 

Recursion formulae and the example of application. After substitution (E4l) - (l25l) into 
when (i y£ k), we have: 



m 



Vi - E 1 6 i>1 



fc-l i 

-k + i)ct l + Y.H c i a . 

j=l p=0 



k-j 
i—p 



+2mE k 8 i>k - 1(1 + 1) b~ ifi b~k,2 
In the case % = k, from f)26p and (1271) . we obtain 
E ■ 

E k ■■ 



k > 1. 



(27) 



^ E 1 = Vl 



2N 2 
1 
2m 



l j 

fc-l k 



j=l p=0 



fc > 1, 



(2? 



that, through the the Taylor-series coefficients for the potential function, is 

(L — 3 N 2 ) V 2 



E 
E 3 
E 4 



iV 2 
" 2 m 2 V 2 
N 2 

' 8m 3 V 4 

(3 L (2 - L) - 5 N 2 (5 - 6 L) - 35 iV 4 ) V V4) , 
AT 4 

= o 4 t,5 ((- 5L (2 + 3L) + 7iV 2 (9-2L)+45iV 4 ) l/ 2 l/ 3 + 
(12 - 50 L + 15 L 2 + 35 iV 2 (3 - 2 L) + 63 N 4 ) V V 5 ) , (29) 



E 1 = V 1 , E 2 



(1-3L + 5N 2 )V 3 , 
((3L 2 -5N 2 -7N 4 ) V 2 2 + 



x 



where N = n + I + 1 and L = 1(1 + 1). 

We see that the zero approximation gives the exact solution for the Coulomb problem. 

As an example of application, we consider the case of the Debye potential, which is 
widely used in many branches of physics: 



a 



V(r) = exp(— fix). 

For this potential, the first corrections to the energy eigenvalues take the form 



(30) 



En 



Ex 



Ea 



E, = 



ma 

N 2 

12 m 2 a 
N 2 

192 m 3 a 2 
N 4 



E l 



a k 



Eo 



(L-3N 2 ) 2 

A K 

Am 



(1 -3L + 5N 2 ) 



K 



(3L (2 + 5L) -5 (11 -6L) N 2 + 77 N 4 ) k 4 , 
(4 - 50 L - 45 L 2 + 35 (7 - 2 L) N 2 + 171 A^ 4 ) 



fi- 



(31) 



320 m 4 a 3 

where N = n + I + 1, and L = 1(1 + 1). 

And again we recognize the explicit perturbation expansion in powers of the small 
parameter k. 

Typical results of the calculation with these formulae are presented in the Table I 
where the sequences of the partial sums of K corrections to the energy eigenvalues for 
the Debye potential V(r) = — aexp(— fix)/r is compared with results of the numerical 
integration, E QVLm , in Coulomb units h — m — a — 1. It is seen that for small values 
of the screening parameter, the convergence of the series is quite sufficient for the use 
them without any renormalization. However, the results become gradually worse when 
the screening parameter increase, as it was pointed for such potentials in 



K 


n = 0, I = 0, k = 0.2 


n = 1, I = 0, k = 0.04 


n = 1, I = 1, « = 0.02 





0.5000000000 


0.1250000000 


0.05555555556 


1 


0.3000000000 


0.0750000000 


0.03555555556 


2 


0.3300000000 


0.0825000000 


0.03805555556 


3 


0.3260000000 


0.0816250000 


0.03781555556 


4 


0.3271000000 


0.0818140625 


0.03786145556 


5 


0.3266800000 


0.0817559375 


0.03784969436 


10 


0.3268179839 


0.0817715528 


0.03785241171 


15 


0.3268059572 


0.0817711705 


0.03785238868 


20 


0.3268100537 


0.0817711991 


0.03785238922 


25 


0.3268067333 


0.0817711951 


0.03785238920 


E 

-^imm 


0.3268085112 


0.0817711958 


0.03785238920 



5 Summary 

In conclusion, a new useful technique for deriving results of the logarithmic perturbation 
theory has been developed. Based upon the fi-expansions and suitable quantization con- 
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editions, new handy recursion relations for solving the bound-state problem for a spherical 
anharmonic oscillator and a static screened-Coulomb potential have been obtained. These 
relations can be applied to excited states exactly in the same manner as to ground states 
providing, in principle, the calculation of the perturbation corrections of large orders in 
the analytic or numerical form. Besides this remarkable advantage over the standard 
approach to the logarithmic perturbation theory, our method does not imply knowledge 
of the exact solution for the zero approximation, which is obtained automatically. And 
at last, the recursion formulae at hand, having the same simple form both for the ground 
state and excited states, can be easily adapted to applying any renormalization scheme 
for improving the convergence of obtained series, as it is described in [32]. 

This research was supported by a grant N 0106U000782 from the Ministry of Education 
and Science of Ukraine which is gratefully acknowledged. 
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